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5 ( Abstract 

We propose to compute the action and global charges of the asymptotically de Sitter solutions in 
Einstein-Gauss-Bonnet theory by using the counterterm method in conjunction with the quasilocal for- 
q ' malism. The general expression of the counterterms and the boundary stress tensor is presented for 

spacetimes of dimension d < 7. We apply this tehnique for several different solutions in Einstein-Gauss- 
Bonnet theory with a positive cosmological constant. Apart from known solutions, we consider also d — 5 
^q' vacuum rotating black holes with equal magnitude angular momenta. These solutions are constructed 

i i , numerically within a nonperturbative approach, by directly solving the Einstein-Gauss-Bonnet equations 

with suitable boundary conditions. 

(N ■ 

> : 

1 Introduction 

On . 

. One of the most fruitful approaches in computing conserved quantities in general relativity is to employ the 

quasilocal formalism pQ. The basic idea here is to enclose a given region of spacetime with some surface, 
and to compute all relevant (conserved and/or thermodynamic) quantities with respect to that surface. For 
qq ' a spacetime that is either asymptotically anti-de Sitter (AdS) [2J, [3], [I] or asymptotically fiat [5], [5J, [7] 

it is possible to extend the quasilocal surface to spatial infinity without difficulty, provided one incorporates 
appropriate boundary terms in the action to remove divergences. The boundary terms are built up with 
• i-H , curvature invariants of the boundary metric and thus obviously they do not alter the bulk equations of 

motion. Therefore, this approach has the nice feature that it is not necessary to embed the boundary 
geometry in a reference background (see also for other applications of this formalism for a different 

asymptotic structure of spacetime). 

The situation is much more involved for asymptotically de Sitter (dS) spacetimes, because of the absence 
of spatial infinity and a globally timclikc Killing vector in this case. In the prescription proposed in 8|, these 
obstacles are avoided by computing the quasilocal tensor of Brown and York pQ (augmented by suitable 
boundary counterterms), on the Euclidean surfaces at future/past timelike infinity I^. This allows also a 
discussion of the thermodynamics of the asymptotically dS solutions outside the cosmological horizon, the 
boundary counterterms regularising the (tree-level) gravitation action as well. The efficiency of this approach 
has been demonstrated in a broad range of examples, including configurations with gravitating matter fields 

cum 

The results in [Sj-[H] concern the case of Einstein gravity with a positive cosmological constant. However, 
for a spacetime dimension d > 4, the Einstein gravity presents a natural generalisation - the so called 
Lovelock theory, constructed from vielbein, the spin connection and their exterior derivatives without using 
the Hodge dual, such that the field equations are second order [16], [17]. Following the Ricci scalar, the 
next order term in the Lovelock hierarcy is the Gauss-Bonnet (GB) one, which contains quadratic powers 
of the curvature. As discussed in the literature, this term appears as the first curvature stringy correction 



1 



to general relativity [181 119j , when assuming that the tension of a string is large as compared to the energy 
scale of other variables. 

In principle, there are no obstacles in computing the action and global charges of EGB solutions in dS 
spacetime by using a quasilocal formalism similar to that proposed in [8] for the Einstein gravity. At any 
given dimension one can write down only a finite number of counterterms that do not vanish at future/past 
timelike infinity. This feature does not depend upon the bulk theory is Einstein or GB. However, the 
presence in this case of a new length scale (the GB coupling constant) implies a complicated expression for 
the coefficients of the boundary counterterms and makes the procedure technically much more involved. 

The corresponding problem for an asymptotically AdS spacetime has been discussed in the recent paper 
[20] (see also [21]). The main purpose of this work is to generalize the boundary counterterms and the 
quasilocal stress energy tensor there to a positive value of the cosmological constant and thus to extend 
the prescription in [8], [12] to the case of EGB theory Our results are valid for configurations with d < 7, 
although a general counterterm expression is also conjectured. In the second part of this paper we apply this 
general formalism to several different asymptotically dS black holes. Apart from known static solutions, we 
consider also rotating black holes with two equal magnitude angular momenta in EGB theory formulated in 
five spacetime dimensions. These solutions are constructed numerically within a nonperturbative approach, 
by directly solving the EGB equations with suitable boundary conditions. They also provide a nontrivial 
generalization in EGB theory of a particular class of the known Myers-Perry-dSs black holes 22J. 

Most of the notation and sign conventions used in this paper are similar to those in ref. [T2] . 



2 The general formalism 
2.1 The action and field equations 

We consider the EGB model with a positive cosmological constant A = (d — 2)(d — l)/2£ 2 , coupled with 
some matter fields with a lagrangean density Cm 

I= TT7^f d d x^9 ( R- 2A+ j L GB + Cm), (2.1) 
16nG J M V 4 / 

where 

Lgb = R 2 — / kR iiV R iXV + RnvarR^^ 7 (2-2) 

is the GB term, while R, R^ u and R^a-r are the Ricci scalar, the Ricci tensor and the Riemann tensor 
associated with the bulk metric g^ v . For d — 4, Lqb is a topological invariant and thus does not contribute 
to the equations of motion; in higher dimensions it is the most general quadratic expression which preserves 
the property that the equations of motion involve only second order derivatives of the metric. The constant 
a in (j2.ip is the GB coefficient with dimension (length) 2 and is positive in the string theory. We shall 
therefore restrict in this work to the case a > 0, although the counterterm expression does not depend on 
this choice. 

The variation of the action (|2.ip with respect to the metric tensor results in the gravity equations of the 
model 

1 a 

Rfiu — -^Rdiiv + Ag^^ + -^H^u = 2T^ V , (2.3) 

where 

H^v = 2(i? MCTKr -R„°^ T — 2R fipua R pa — 2R tia R° ' v + RR^) — -£gb3^ > (2-4) 

and T^ u is the energy-momentum tensor of the matter fields. 

For a well-defined variational principle, one has to supplement the action (|2.ip with the Gibbons-Hawking 
surface term [33] 

li E) = - o-p / d d - l x^jK , (2.5) 
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and its counterpart for the GB gravity 



^ = - / d d -'x^ (sJ - 2G ab K ab ) , (2.6) 



where 7 a h is the induced metric on the boundary with the outward-pointing normal vector n° and K is the 
trace of the extrinsic curvature K ab of the boundary. The factor s in (|2.6[) is s = +1 for a spacelike normal 
vector n a and s = — 1 for a timelikc normal vector (the case considered here), see e.g. |24j . Other quantities 
in (|2.6[) are G ab - the Einstein tensor of the metric j a b and J -the trace of the tensor 

J ab = X -{2KK ac K c h + K cd K cd K ab - 2K ac K cd K db - K 2 K ab ) . (2.7) 

Also, the case of interest in this paper corresponds to a spatial boundary future/past timelike infinitjQ. 

Therefore dAi^ are spatial Euclidean boundaries at early and late times, while Jg^- d d ~ 1 x indicates an 
integral over the late time boundary minus an integral over the early time boundary. In what follows, to 
simplify the picture, we will consider the I + boundary only, dropping the ± indices (similar results hold for 

The equations (|2.3[) present many interesting solutions possessing new features as compared to the pure 
Einstein gravity case (for a review, see the recent work [27] ). In this Section we shall consider the issue of 
computing the action and global charges of the EGB solutions in dS spacetime, several examples of such 
configurations being discussed in the next Sections. Here we have found convenient to wr 



-j—, with U=^ f L + l, (2.8) 

which results in a compact form for the expression below; physically, £ c corresponds to an effective dS length 
scale in EGB theory. 

2.2 The counterterms and the boundary stress tensor 

The GB term in (|2.ip does not change the general formalism to compute the conserved charges and the 
action of asymptotically dS solutions developed in [8], [12]. Therefore, we only recapitulate the basic steps 
here, emphasizing the new features which emerge for a^O. 

In general, the action (|2.1[) (together with the boundary terms (|2.5[) . (|2.6|) ) is divergent when evaluated 
on a solution to the equations of motion (|2.3j) . In the counterterm approach, the remedy is to supplement 
the initial action (12. 1|) by a boundary counterterm part I c t depending only on geometric invariants of the 
boundary metric (therefore the bulk equations of motion remain the same). I c t regularizes the tree gravita- 
tional action I c \ and the boundary stress tensor. Crucial to the success of the counterterm prescription is 
that the divergencies are universal, so that a single choice of the counterterms suffices to render finite the 
action of all asymptotically dS solutions. 

For d < 8 solutions, our proposal for the boundary counterterm action itd 

1 ' <-w- ( ^)(^) + ^m?p-ioR M 



8ttG J dM v '{ £ c 3 ' 2(d-3) 
£ 3 c Q(d-6) 

~ 2(d - 3) 2 (d - 5) f V " 4 ( d - 2 ) " J 2 ( d - 4) 



u( Ra 6 R Qb - 77^4^ - Jt^(U - 1)Lgb 



1 For the black hole solutions in this paper, this corresponds to evaluate various quantities for some fixed radius larger than 
the radius of the cosmological horizon and then sending this radius to infinity. 

2 For the sake of simplicity, we have restricted ourselves to the case of solutions with a well defined Einstein gravity limit. 
However, the results in this section can easily be generalized to the branch of solutions diverging as a — ► 0, in which case several 
terms in l|2.9|l have an opposite sign. 

3 Note that in odd spacetime dimensions, for some boundary geometries, there is an additional logarithmic divergence that 
cannot be cancelled without includding an explicit cutoff dependence in the counterterm action, which should be supplemented 
with new extra terms. This feature occurs already for Einstein gravity theory, leading to a conformal anomaly similar to what 
has been obtained in the context of AdS spacetime [2]. However, this is not the case of the solutions discussed in next sections. 
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where R, R 06 and Lqb are the curvature, the Ricci tensor and the GB term associated with the induced 
metric 7. Also, Q(x) is the step-function with (x) = 1 provided x > 0, and zero otherwise. One can 
easily see that as a — > (thus U — > 1), the known counterterm expression in the Einstein gravity [5], |12j is 
recovered. 

Following ref. [20] . we conjecture the general form of I ct in d spacetime dimensions: 




^ V7 E 6 ( d - 2fc ) h( u ) L E + f2{U)C(k-i) \, (2-10) 



where Le is the corresponding fc-th part of the counterterm lagrangean for a theory with only Einstein 
gravity in the bulk (with the length scale I in front of it replaced by the new effective dS radius £ c ) and 
£(fc_i) is the (fc — 1) term in the Lovelock hierarchy. The functions fi{U), f2(U) are first order polynomials 
in U, whose expression can easily be derived from those given in [5D] for A < 0. The series (|2.10[) truncates 
for any fixed dimension, with new terms entering at every new even value of d. 

Once we know the expression of the boundary counterterm, the computation of the conserved charges 
is performed in a similar way to the a = limit [H [12] . The (Euclidean) boundary metric on equal time 
surfaces can be written, at least locally, in a ADM-like general form 

ds 2 = lab dx a dx b = N 2 p dp 2 + a ab (diP a + N a dp) (dip b + N b dp) , (2.11) 

where N p and N a are the lapse function and the shift vector respectively, while ip a are angular variables 
parametrizing a closed surfaces E. The physical significance of the coordinate p in (|2.11[) depends on the 
considered situation; e.g. for the black hole solutions discussed in the next Sections, p is the coordinate 
associated with the asymptotic Killing vector that is timelike inside the static patch of dS, but spacelike at 
1±. 

Varying the total action with respect to the boundary metric jab results in the following boundary 
stress-energy tensor 

j, b= ±^(, + ,r + ,r> + , a ), (2.i2) 

with the following expression valid for d < 8: 
87rGT ab = K ab - lab K + |(Q o6 - l -Q lab ) + ^ 7ab (^±^) + ^£^_il( 2 - U) (^R ab - l - lab R 

+e " Q ^ ~ 6 ) { " ( d -3)V5) ( " ^ -W^) R2 )- W^) RRab + 2RCdRcadb 

^W^2)^ WbR + y2R - ~ w-i)^ v2r ) + 2 { d-m-\ ){ d-,M + • ■ - (2 - 13) 

where [H], [H] 

Q ab - s (2KK ac K r b - 2K ac K cd K db + K ab {K cd K cd - K 2 )) + 2KR ab + RK ab - 2K cd R cadb - 4R ac K b c , (2.14) 

and H ab given by (|2.4[) , this time written in terms of the boundary metric j abl however. All terms in (|2 . 13[) . 
except the first four, come from the variation of the counterterms in (|2.9j) . The boundary stress-energy 
tensor T a b measure the response of the spacetime to changes of the boundary metric and encodes the notion 
of conserved global charges. 

Following [5], [T2], let us suppose that is a Killing vector generating an isometry of the boundary 
geometry (|2. 1 1|) . Then it is straightforward to show that T^- 7 is divergenceless and one can define a 
conserved quantity 0£ associated with £ l as follows 

£5 = j <l"r s n „■!,::■. (2.15) 
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where n l is a unit vector normal on a surface of fixed p. The physical interpretation of this relation is the 
same for any theory of gravity: it means that a collection of observers, on the hypersurface with the induced 
metric cr^, would all measure the same value of 0£ provided this surface has an isometry generated by £ l . 
As mentioned above, a dS spacetime has no globally timelike Killing vector, which makes difficult to define a 
mass for the solutions with this asymptotics. However, for all cases of interest {e.g. the black holes solutions 
in the next Sections), there is a Killing vector that is timelike inside a static patch, while it is spacelike 
outside the cosmological horizon and therefore at T ± . (Moreover, any spacetime that is asymptotically dS 
will have such an asymptotic symmetry generator.) The total mass/energy of solutions is evaluated with 
respect to this Killing vector. 

Proceeding further, one can define a Hawking temperature Tu (e.g. by computing the corresponding 
surface gravity) and entropy S for the cosmological horizon by using the saddle point approximation to the 
gravitational partition function (namely the generating functional analytically continued to the Euclidean 
spacetime). In the semiclassical approximation, the dominant contribution to the path integral will come from 
the neighborhood of saddle points of the action, that is, of classical solution; the zeroth order contribution 
to logZ is given by —I c i- A tree-level evaluation of the path integral with a GB term may be carried out 
along the lines described e.g. in ref. |12j for the Einstein gravity case. Therefore, we find the entropy of the 
cosmological horizon (with [3 = 1/Tjj) 

S = P(E - md) - Id, (2.16) 

which is found upon application of the Gibbs-Duhcm relation to the partition function, with chemical 
potentials £j and conserved charges /Uj, while E is the total mass/energy, evaluated according to (|2.I5I) . 
Also, all solutions should satisfy the first law of thermodynamics for the cosmological horizon 

dS = P(dE - mdU), (2.17) 

which provides a test of the general formalism. 



3 Applications: known solutions 

3.1 dS spacetime in EGB theory 

As the simplest illustration of the above formalism, we consider the case of empty dS spacetime. This solution 
has a simple form in a large number of coordinate systems. For example, there is a static frame centered 
on each observer (timelike geodesic) in dS. Moreover, when a black hole exists, there is still a static frame 
centered about the black hole. Since different parametrizations emphasize different features, it is of interest 
to consider dS spacetime in alternative coordinate systems corresponding to different classes of observers. 
Starting with an inflationary coordinate system, the dS solution reads 

ds 2 = -dt 2 + e 2t '^dx 2 , (3.1) 

which solves the EGB equations (|2.3|) with T M „ = (i.e. no matter fields). The properties of this solution 
are similar to the case of Einstein gravity (see e.g. [26)): the equal time surfaces here are flat, while t runs 
from — oo to +oo. One can easily verify that the counterterms (|2.9p removes all divergencies of the total 
action for d < 7, and leads to I c i = 0. The total mass/energy of this solution is also vanishing, since Tjj = 0. 
The situation is different for a static coordinate system, the corresponding line element being 

dr 2 

ds 2 = —-+r 2 dn 2 _ 2 ~F(r)dt 2 , (3.2) 
F(r) 

where (here we shall consider only the branch of solutions with a smooth Einstein gravity limit) 



F{r) = 1 + a( d -3)(d-4) V - V 1 + -^ d -^)- I 2 
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This spacetime has a cosmological horizon at r c = £ c (where F(r c ) = 0), with an associated temperature 
Tfj = (2tt£ c )^ 1 . The topology of this solution for large constant r > r c , is an Euclidean cylinder R x S d ~ 2 
and t is the coordinate along the cylinder. X ± are located outside the future/past cosmological horizons, 
where r is timelike and t is spacelike. The relationship between the coordinate patches (13.1[) and (|3.2[) and 
their Penrose diagrams are presented in ref. [55] . 

The general formalism in Section 2 is applied working outside of the cosmological horizon, where F(r) < 0. 
The gravitational mass/energy is the charge associated with the Killing vector d/dt — now spacelike outside 
the cosmological horizon. As expected, the total energy found by using the counterterm prescription vanishes 
for an even dimensional spacetime and has a nonzero value for an odd d: 

V d - 2 (d-2)!! 2 f(d-2)U-2\. d _ 3 > 



Mo = 8^7^1)7 ^ d-4 ) ****** . (3-3) 

where Vd-i is the area of the unit S d ~ 2 sphere and p > 2 is an integer. For solutions in Einstein gravity 
(a = 0), Mo is usually interpreted as the Casimir energy in the context of dS/CFT correspondence. Also, it 
reduces to the expression obtained in ref. p~2j when [7 = 1. 

From (|2.16p one finds the following expression for the entropy of dS spacetime in EGB theory: 



S = 



V d . 



e 4 (^c + f(rf-2)(d-3)), (3.4) 



4G 

which in the limit of small a can written in the simple form 

S = S + S c with So = £ ,^~ 2 , S c = c^-^d{d - 2)(d - 3)£ d ~ 4 . (3.5) 

From the study of (|3.1|) . ()3.2|1 we conclude that, similar to the case of Einstein gravity, the horizon and 
entropy of the dS space in EGB theory have an obvious observer dependence. 



3.2 Reissner— Nordstrom-dS-GB black hole 

These solutions are found for a matter lagrangean density Cm = —J~ 2 , with the Maxwell field strength 
tensor T = dA, where the (pure electric-) gauge potential is 

where $ is a constant. Working again in a static coordinate system, the line element is still given by (|3.2[) , 
with a different expression for F(r), however: 

2r 2 / / M CP T~ \ 

^ - 1 + a( d-3)(d-4) V - V 1 + ^ ^ 3)(d ~ 4)( ^ - ^) + P ] ) ' (3J) 

As argued below, M and Q in the above expression are constants proportional to the gravitational mass /en- 
ergy E and the total electric charge Q, respectively. The Q = limit of this metric corresponds to the EGB 
generalization of the McVittie solution describing a Schwarzschild black hole embedded in dS spacetime [2] . 

A discussion of the solution (|3.6[) , (|3.7|) appeared in ref. [29] (see also ref. [13] for an extended analysis 
of the limiting case a — 0, including also multi-black hole generalizations). Here we briefly review its basic 
properties. One can easily verify that the metric has a curvature singularity at the origin r = 0. In general, 
the metric p.2[) presents Killing horizons at the radii where F(r) vanishes. Of interest are the outer black 
hole horizon at r = r/, and the cosmological horizon r — r c corresponding to the largest root of F(r). The 
Hawking temperature of associated to each of the horizons is To* = |F'(r/ l c )|/(47r), where a prime denotes 
the derivative with respect the radial coordinate. The two horizons are not in thermal equilibrium because 
the time periods in the Euclidean section required to avoid a conical singularity at both do not match in 
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general. An extremal black hole is found by imposing F{rh) — F'(i"h) = which fixes M, Q as functions of 
t, a and r/j (a similar relation is found when considering instead the cosmological horizon) . The constant 
in (|3.6[) is usually fixed such that A t (r c ) = 0, and thus it corresponds to the electrostatic difference between 
the cosmological horizon and infinity. 

The computation of the mass, action and entropy of a RNdS black hole is a direct application of the 
method described in the previous section. The gravitational mass/energy is the charge associated with the 
Killing vector d/dt. The total mass/energy found by using the counterterm prescription described in the 
previous Section is 

E = -^(d-2)M + M , (3.8) 

with M the Casimir term given by (|3.3p . The negative sign implies that the black hole lowers the total bulk 
energy with respect to the total energy of the pure dS spacetime [8]. 

The computation of the total electric charge is similar to that performed in [13] for the case without a 
GB term. The results there show that the total electric charge evaluated at future/past infinity is 

Q = ^V2(d-3)(d-2). (3.9) 

From (|2.16p one finds the entropy of the comological horizon (note that both So and S c have a nontrivial 
dependence on a): 

S = S + S C with S a = ^rt 2 , S c = a^K d -\d-2)(d-Z). (3.10) 

One can easily verify that the first law of thermodynamics (|2.17[) also holds, with /ij = <f>, £j = Q. 

It would be interesting to study the properties of the Reissner-Nordstrom-de Sitter solution (|3.6|) , (|3.7p 
in the inflationary coordinate system (|3.ip . This problem has been considered in ref. [13] in the absence 
of a GB term in the action. Interestingly, the same general picture have been found there for black holes 
in both coordinate systems (|3.ip and ()3.2|) . which shows the complex relation between different classes of 
observers in dS spacetime. For example, the mass of the black holes in an inflationary coordinate system is 
still given by f3~% with M = howeveiQ. We expect that a similar result will be found in the presence 
of a GB term. However, in the absence of an explicit form of the Reissner-Nordstrom-dS-GB black hole in 
the inflationary coordinate systerr[f|, any progress in this direction appears to require a separate numerical 
study of these solutions. 



4 Rotating EGB black holes with positive cosmological constant 
4.1 The metric ansatz and known limits 

The computation of the global charges and entropy of a rotating black holes in EGB theory represents 
another nontrivial aplication of the general formalism in Section 2. Unfortunately, no exact solutions are 
available in this case, and one has to solve numerically the field equations. 

To simplify the general picture, we consider here the vacuum case in d = 5 dimensions only, although the 
inclusion of a U(l) field is straightforward in principle. A general spinning black hole solution is characterized 
in this case by two angular momenta and its mass/energy, and can be found by solving a set of seven partial 
differential equations. However, the numerical problem is greatly simplified by taking the apriori independent 
two angular momenta to be equal in order to factorize the angular dependence [31] . [32] . The asymptotic 
expressions and the explicit computation of the action and boundary stress tensor also simplifies drastically 
in this case. 

4 This result has been interpreted in 1131 as providing support for the putative dS/CFT correspondence, since the general 
features of the CFT dual to a black hole should not depend on the dS slicing choice. 

5 The main obstacle here is the absence of a simple closed form expression of the Reissner-Nordstrom-GB (or even 
Schwarzschild-GB) black hole in an isotropic coordinate system for the A = case. For a = 0, this form of the solution 
is used to construct cosmological configurations by using the prescription in [301 - 
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To construct these solutions, we use the same metric ansatz employed in ref. [33j for the corresponding 
problem with A < 0: 

ds 2 = 4^7 + g(r)d6 2 + h(r) sin 2 6 (dp - w{r)dt) 2 + h(r) cos 2 9 (dib - w(r)dt) 2 (4.1) 
fl r ) 

+ (g(r) - h(r)) sin 2 9 cos 2 0(dip - d^) 2 - b(r)dt 2 , 

where 9 £ [0, 7r/2], (p, -0) £ [0, 2tt], r and t being the radial and time coordinates. This ansatz has a residual 
degree of freedom which is fixed by taking g(r) = r 2 . 

The equations satisfied by the functions b,f,h,w result directly from (|2.3p . We refrain to write them 
because they are very long and not particularly enlightening. They present however two exact solutions 
which are important in what follows. The dSs generalization 22\ of the Myers-Perry rotating black holes 
[51] with equal magnitude angular momenta (hereafter MPdSs) is found for a — (no GB term) and has 

t 2 ZMZ 2Ma 2 , / 2Md 2 \ , , 2Ma .. 2 ,.. r 2 f(r) /A . 

f( r ) = — + — > Hr) = r 2 (l + —J , w(r) = g{r) = r 2 , 6(r) - (4.2) 

where M and a are two constants related to the solution's mass and angular momentum, while 5 = 1 + a 2 /£ 2 . 

For g(r) = h(r) = r 2 , w(r) = and f(r) = b(r) = 1 + r 2 /a^l - + 2a(M/r i + l/£ 2 )^j , one 

recovers the Schwarzschild-dSs solution with a Gauss-Bonnet term. The slowly rotating generalisation of 
this solution^! is found for small values of the rotation parameter a, and reads 



2aU 2 / / 2MR(U -\) \ fl n . 

W{r) = mr^) V 1+ w - 1 ' (4 - 3) 



the other metric function remaining unchanged to this order in a. 



4.2 Boundary conditions and global charges 

We want the generic line element (|4.1[) to describe a nonsingular, asymptotically de Sitter spacetime outside 
a cosmological horizon located at r = r c > 0, with f(r c ) = 0. Here f(r c ) = is only a coordinate 
singularity. The regularity assumption implies that all curvature invariants at r — r c are finite. Outside 
the cosmological horizon r and t changes the character {i.e. r becomes a timelike coordinate for r > r c ). 
A nonsingular extension across this null surface can be found just as at the event horizon of a black hole. 
These configurations possess also an event horizon located at some intermediate value of the radial coordinate 
< Th < r c , all curvature invariants being also finite as r — > ru- 

Restricting to nonextremal solutions, the following expansion holds near the event horizon with the 
parameters f^, b^, w\ and h%, where (/f - , 6^, h%) > 0: 

f(r) = fHr-r h )+0(r-r h ) 2 , h(r) = h h h + 0(r - r h ), (4.4) 
b(r) = b1(r - r h ) + 0(r - r h ) 2 , «,(r) = w\\ + 0(r - r h ). 



A similar expansion holds for cosmological horizon, the corresponding parameters there being /°, bl, h' 



and Wfr. 



Both the event and the cosmological horizon have their own surface gravity n h,c , the associated Hawking 
temperatures being 



■A,c \^ H " 



H 2?r 4tt v ' 



6 By using the results derived in |35| for A < 0, one can find a different set of closed form EGB asymptotically dS^ rotating 
black hole solutions with only one nonvanishing angular momentum (where the rotation parameter appears as a small quantity), 
the effects of an U(l) field being also included. A different approximation of the rotating black hole solution of the d = 5 EGB 
equations with a cosmological constant have been presented in closed form in [36] ■ 
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Another quantities of interest are the area A^ c of the black hole/cosmological horizon 

A h H c = Vh~ h h ' C (r h h c ) 2 V 3 , (4.6) 

where V3 = 2ir 2 denotes the area of the unit three dimensional sphere. 

The Killing vector x — d/dt + fl v d/dip + VL^d/dip is orthogonal to and null on both horizons. For the 
solutions within the ansatz (|4.1|) . the event horizon's angular velocities are all equal, f^' c = f^ ,c = £l^ c = 
w(r)\ r = rhc . 

A direct computation reveals that the solution admits at large r a power series expansion of the form: 

/«=i-s+!>^r. *)-.-£ + g-(£)" <4j) 

2fc \ -, / n \ 1k+l 



_ W 2 k+1 — 
fe>l x k>l v 

where the coefficients f 2 k,b 2 k, h 2 k,w 2 k+i with k > 1 are determined by f 2 ,b 2 and W3. Specifically, we find 
U = h i = 62-/2, &6 = (6 2 (/ 2 + 6 2 (i7-2)) + (3[/-2)w|)/(2i7), w 7 = -(/a + U(2b 2 - 3f 2 ))w a /(2U), for the 
lowest order nonvanishing terms. Their expression becomes more complicated for higher k, with no general 
pattern becoming apparent. 

The mass/energjjj E and angular momenta of these solutions evaluated at future/past timelike infinity 
by using the counterterm formalism are fixed by the constants f 2l b 2 and W3, and read 



The entropy of these solutions associated with the cosmological horizon is found from the relation (|2.17|) 
with m = Q,^, € = J: 



S = S Q + S GB , with ^0 = 7^, S GB = a-±^ h (4:-—^). (4.9) 



4.3 The numerical method 

Finding numerical solutions of a field theory model in a dS spacetime is a notoriously difficult task. Therefore, 
before describing the properties of the solutions, we shall give some details on the numerical methods we have 
usecj^). The EGB field equations were solve by employing a collocation method for boundary-value ordinary 
differential equations, equipped with an adaptive mesh selection procedure [37]- Typical mesh sizes include 
fO 3 — fO 4 points. The solutions have a typical relative accuracy of 10~ 8 . In constructing rotating EGB-dS 
black holes, we make use of the existence of the MPdSs and Schwarzschild-GB-dS closed form solutions, and 
employ them as starting configurations, increasing gradually f2^ or a, respectively. 

However, when trying to find black hole solutions with A > for r G [r^oo] by imposing a regular 
horizon at r — r/i, one has to tackle the technical difficulty that there also appear a cosmological horizon^. 
That is, the metric functions /, b admit a zero at an intermediate value of the variable r, say at r = r c > rh- 
Of course, the value of r c is not known apriori as a function of A, a. However, the condition of a regular 
horizon should be imposed both at r — and r = r c . In our approach, we impose by hand the values 
of r/j,r c and solve the equations first for r € [r/,,r c ] as a boundary value problem. At the same time, we 
compute the value of A corresponding to this cosmological horizon by using the fictious equation dk/dr = 0. 

7 In the expression of E, we have subtracted the Casimir energy of the pure dSs space as given by 1 13. 31 . 

8 The our knowledge, this is the first attempt in the literature to numerically construct EGB rotating solutions in a dS 
background. The approach and the numerical methods here are quite different from those employed e.g. in I33| for rotating 
EGB solutions with AdS asymptotics or for rotating Einstein-Maxwell black hole solutions 1311 . 1321 . 

9 We do not consider in this work the behaviour of solutions inside the black hole even horizon r < r^. 
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In a second step, we finally integrate the equations for re [r c , oo] as an initial value problem with this value 
of the cosmological constant. This assures that the metric functions and their first and second derivatives 
are continuous at r = r c . 

In this approach, the set of boundary condition we imposed at rh,r c is 



/ = 0, 6 = 0, b' = 1, G{gij,g' tj ) = 0, w = Wh at r = r h , 
and f = , 6 = 0, G(gij , g'^) = at r = r c . 



(4.10) 
(4.11) 



where G(gij, g[j) is a complicated expression in terms of the metric function and their first derivatives which 
occurs from the condition for a regular horizon. In the above expression, the arbitrary rescaling of time is 
used to set b'(rh) — 1, keeping in mind that the function 6(r), w(r) have to be renormalizcd at the end of 
the second step according to 



6(r) — > b(r) = 6(r)/i 2 , w(r) —> w(r) = w(r)fj,, 



(4.12) 



where the constant /i is chosen in such a way that the 6(r) approaches the asymptotic (|4.7[) . 

One disadvantage of this method is that the solutions cannot be studied systematically for fixed A. For 
the same reason, we have found difficult to study families of solutions obtained by varying a while fl H is 
fixed. 

To summarise, in our approach the input parameters are the black hole event horizon radius rh, the 
cosmological horizon radius r c , the black hole event horizon velocity f2^ and the GB coupling parameter a. 
The value of the cosmological constant, the metric functions and their derivatives at r = r c and the global 
charges emerge from the numerical output. 




0.4 

logioto 

Figure 1: The profiles for a generic d — 5 rotating black hole solution in EGB theory with positive 
cosmological constant. 



4.4 Numerical results 

A systematic study of the properties of these rotating black holes appears to be a difficult task and is beyond 
the purposes of this work. In practice we have solved the equations numerically for several values of r^, r c 
and flff and a range of the Gauss-Bonnet coupling constant a. 
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Figure 2: Several parameters are compared for the Myers-Perry-dSs (dotted curves) and Einstcin-Gauss- 
Bonnct (continuous line) rotating black hole solutions. The input parameters here were r^, = 1, r c = 3, 
a = 0.5 and t 2 = 10. 

When increasing from zero the angular velocity Q H , we have found numerical evidence for the existence 
of nontrivial generalizations of any static Schwarzschild-GBdS configuration we considered; the shape of 
the metric function w(r) we found for small values of VI H is in good agreement with (|4.3[) . We reach the 
same conclusion when considering instead GB counterparts of the Einstein gravity rotating solution (|4.2p . 
by slowly increasing the parameter a. As a general remark, the qualitative features of all solutions we have 
constructed are rather similar to the MPdSs case. For a > 0, we have noticed only quantitative difference 
in the values on the cosmological horizon and at infinity, for a number of parameters of interest. 

In order to limit the amount of numerical investigation, we have studied in details mainly the case 
r/j = 1, r c = 3. For the non-rotating limit, this corresponds to i 2 = 10 and — (&ttG)E /Vj, = 3M/4 = 
3(5a + 9)/40, 4irT H = 8/(5(a + 1)). 

The profiles of the metric functions of a typical EGB-dS rotating black hole solution corresponding to 
a = 1, il H = 0.66, i ~ 3.2 are presented on Figure 1. One can see that the rotation leads to non constant 
values for h(r)/r 2 and b(r) ^ f(r), especially in the region close to the black hole horizon. Also, the metric 
functions and their derivatives are continuous at the cosmological horizon (although to simplify the plot we 
presented there only the profile of b'{r)). 

Several parameters characterizing the solutions are represented in Figure 2 as a function of the angular 
velocity at the black hole horizon. The data corresponding to MPdSs solution is represented by the dotted 
lines and results from analytical calculations. In contrast, the curves correponding to the EGB theory are 
represented by continuous lines and result from our numerical calculation. (The energy E and angular 
momentum J are represented in the units of V^/{At:). The Hawking temperature is represented in units 
l/(47r) and horizon area Ah in units V3, while we have set also G = 1 in all data.) Along with the case of 
the MPdSs solutions, the EGB black holes exist up to a maximal value of Sl H = £l m ax- For a = one finds 



For the cases we have investigated, when fixing the values of r^, r c , this maximal value gets larger when 
the Gauss-Bonnet coupling constant increases. In this limit, the solution approaches an extremal black hole, 




(4.13) 
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Figure 3: The black hole event horizon area is plotted as function of the black hole temperature for several 
values of the GB coupling constant a. These solutions have = 1, r c = 3 and £ 2 = 10. 

i.e. the functions f(r),b(r) reach a double root at the black hole event horizon. However, these function still 
present a simple zero at the cosmological horizon. For the values of the parameters adopted here we find 
e.g. Q h H « 0.64 for a = and fl 1 ^ w 0.72 for a = 1. 

For the Einstein gravity black holes, the event horizon area increases with increasing while 
the Hawking temperature decreases. The entropy S(Th) turns out to be a decreasing function of the 
temperature, for fixed event horizon radius. We have found that when the GB parameter a is large enough, 
the scenario is completely different. For instance, both Tk and A\ decreases while fi^ increases. The 
entropy is an increasing function of T^. In Figure 3 we show this behaviour for several values of a. There 
the parameter varying on the A^(T^) curves is . 

As with other rotating black holes, these solutions present also an ergoregion inside of which the observers 
cannot remain stationary, and will move in the direction of rotation. The ergoregion is the region bounded 
by the black hole event horizon, located at r = rt and the stationary limit surface, or the ergosurface, with 
r = te < r c . The Killing vector d/dt becomes null on the ergosurface, i.e. gtt = —b(rE) + r E w(rE) 2 = 0. For 
the ansatz (|4.1[) . the ergosurface does not interesect the horizon. We observe indeed that, for rotating solution 
with fixed f2^, the value te decreases slightly and get closer to r/, when a increases. For exemple, with 
= 0.66, we get Te/th w 1.34 and r^/r/, ~ 1.27 respectively for a = 0.1 and a = 1. In principle, there 
is also a second value of r, located outside the cosmological horizon, where the Killing vector d/dt becomes 
null again. However, for all solutions we could construct, the metric component g tt there is dominated by 
b{r) and thus the "cosmological" ergo-radius hardly differs from r = r c . 

Another qualitative difference between Einstein and EGB black holes resides in the magnitude of the 
ratio p = f'/b'\ r — rh . For instance, for a = 0, we have p < 1 for all values of the angular momentum. 
However, when a got sufficiently large, one can find solutions with p > 1. 

5 Further remarks 

In this work we have presented the boundary counterterm that removes the divergences of the action and 
conserved quantities of the solutions in EGB theory with a positive cosmological constant for a spacetime 
dimension d < 7. Similar to the case of Einstein gravity, the counterterm is built up with curvature invariants 
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of the boundary metric. Their coefficients, however, present an explicit dependence of the dimensionless 
factor a 2 A. 

Here one should say that the expression of the counterterm proposed in this paper was obtained by 
demanding cancellation of divergencies for a number of solutions in EGB theory, which was also the approach 
used in initial work on the boundary counterterm in Einstein gravity [2J, [3J. However, for asymptotically dS 
solutions in the Einstein gravity, there exist an algorithmic procedure for constructing I ct in a rigurous way, 
and so its determination is unique for a = [12] . This procedure involves solving the Einstein equations 
(written in Gauss-Codacci form) in terms of the extrinsec curvature functional of the boundary and its 
derivatives to isolate the divergent parts. All divergent contributions are independent of the boundary 
normal and can be expressed in terms of intrinsic boundary data. In principle, this approach can be extended 
to asymptotically dS solutions in EGB theory, the only obstacle we can see at this stage being the huge 
complexity of the required computation. A more promising direction would be to look for the expression of 
I c t in the linear order in a, by generalising the work in [21] to the A > case. 

For asymptotically AdS solutions, an alternative regularization prescription for any Lovelock theory 
has been proposed in [38) . This approch uses boundary terms with explicit dependence on the extrinsic 
curvature K a b, also known as Kounterterms. It would be interesting to generalize the approach in [38 to dS 
asymptotics and to compare the results with those found here. 

In the second part of this work, the general formalism has been applied for several different asymptoti- 
cally dS black hole solutions in EGB theory. Apart from several known solutions, we have considered also 
rotating black holes with two equal-magnitude angular momenta in d = 4 + 1 EGB theory with a positive 
cosmological constant. Although the numerical difficulties associated with the existence of a cosmological 
horizon prevented us from a systematic study of the parameter space, we have presented arguments for the 
existence of nontrivial generalization in EGB theory of a particular class of the known MPdSs black holes. 

As avenue for future research, it would be interesting to consider the status of " the maximal mass 
conjecture" in EGB theory, by using the mass definition proposed in this work. Formulated in [5] for 
Einstein gravity, this conjecture states that any asymptotically dS spacetime cannot have a mass larger than 
the pure dS case without inducing a cosmological singularity. Here we mention only the fact that all rotating 
black holes we have constructed in Section 4 satisfy this conjecture. 

The conserved charges of the rotating solutions in this paper have been evaluated on a Euclidean surface 
at future timelike infinity. In principle, by using the the results in Section 2, a similar computation can be 
performed for a spatially finite boundary inside the cosmological event horizon. The corresponding problem 
for Kerr-dS rotating black holes in Einstein gravity has been considered in ref. [39] . The results in that work 
show that quasilocal angular momentum is independent on the radius of the boundary, which does not hold 
for the total mass of the solutions. 

The relevance of the results discussed in this paper in a dS/CFT context is another interesting open 
question. For the AdS/CFT case, the higher derivatives curvature terms can be viewed as the corrections of 
large N expansion of the boundary CFT in the strong coupling limit, see e.g. [40] . For the asymptotically 
dS case, any progress in this direction is likely to require first a better understanding of the conjectured 
dS/CFT correspondence [41] with a = 0. 
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